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Abstract 

We give the distribution functions, the expected values, and the moments of Hnear 
combinations of lattice polynomials from the uniform distribution. Linear combi- 
nations of lattice polynomials, which include weighted sums, linear combinations 
of order statistics, and lattice polynomials, are actually those continuous functions 
that reduce to linear functions on each simplex of the standard triangulation of the 
unit cube. They are mainly used in aggregation theory, combinatorial optimization, 
and game theory, where they are known as discrete Choquet integrals and Lovasz 
extensions. 
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1 Introduction 

Let h : [0, 1]" — M be an aggregation function and let X be a random vector 
uniformly distributed on [0, 1]". An interesting but generally difficult problem 
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is to provide explicit expressions for the distribution function and the moments 
of the aggregated random variable Y = ^(X). 

This problem has been completely solved for certain aggregation functions (see 
for instance pT| §7.2]), especially piecewise linear functions such as weighted 
sums [S] (see also [H]), linear combinations of order statistics pi2Uf25j (see 
also fT'j §6.5] for an overview), and lattice polynomials [TS], which are max-min 
combinations of the variables. 

In this note we solve the case of linear combinations of lattice polynomials, 
which include the three above-mentioned cases. Actually, linear combinations 
of lattice polynomials are exactly those continuous functions that reduce to 
linear functions on each simplex of the standard triangulation of [0,1]"". In 
particular, these functions are completely determined by their values at the 
2" vertices of [0, 1]". 

The concept of linear combination of lattice polynomials is known in combi- 
natorial optimization and game theory as the Lovdsz extension [5|13|15|I2^3J of 
a pseudo-Boolean function (recall that a pseudo-Boolean function is a real- 
valued function of 0-1 variables). When it is nondecreasing in each variable, it 
is known in the area of nonlinear aggregation and integration as the discrete 
Choquet integral [TUII14fl6j . which is an extension of the discrete Lebesgue 
integral (weighted mean) to non-additive measures. The equivalence between 
the Lovasz extension and the Choquet integral is discussed in [T6] . 

This note is set out as follows. In Section 2 we elaborate on the definition 
of linear combinations of lattice polynomials and we show how to concisely 
represent them. In Section 3 we provide formulas for the distribution function 
and the moments of any linear combination of lattice polynomials from the 
uniform distribution. Finally, in Section 4 we provide an application of our 
results to aggregation theory. 

Throughout we will use the notation [n] := {1, . . . Also, for any subset 
y4 C [n], will denote the characteristic vector of A in {0, 1}". Finally, for 
any function h : [0, 1]*^ M, we define the set function Vh '■ 2^"' — M as 
VhiA) := /i(1a) for all A C [n]. 



2 Linear combinations of lattice polynomials 

In the present section we recall the definition of lattice polynomials and we 
show how an arbitrary combination of lattice polynomials can be represented. 

Basically an n-place lattice polynomial p : [0, 1]" ^ [0, 1] is a function defined 
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from any well-formed expression involving n real variables linked 
by the lattice operations A = min and V = max in an arbitrary combination 
of parentheses (see e.g. Birkhoff [6, §11.2]). For instance, 

X2, xs) = {xi A X2) V X3 

is a 3-place lattice polynomial. 

Consider the standard triangulation of [0, 1]" into the canonical simplices 

5^ := {x G [0,1]" I ^ ■■■ ^ ((tG6„), (1) 

where 6„ is the set of all permutations on [n]. Clearly, any linear combination 
of n-place lattice polynomials 

m 
i=l 

is a continuous function whose restriction to any canonical simplex is a linear 
function. According to Singer pSl, §2], h is then the Lovdsz extension of the 
pseudo-Boolean function /?.|{o,i}", that is, the continuous function defined on 
each canonical simplex S^j as the unique linear function that coincides with 
/i|{o,i}" at the n + 1 vertices 

'■= l{<T(i),...,(T(i)} (i = 0, . . . , n) 
of Sa- It can be written as [23, §2] 

n 

h{^) = Y.(h1-hU)x.i^) (xG5.), (2) 

i=l 

where hi := hief^ = Vh(^{o'{l), . . . , cr(i)}^ for alii = 0, . . . , n. In particular, 
K = 0. 

Conversely any continuous function h : [0, 1]" —>■ M that reduces to a linear 
function on each canonical simplex is a linear combination of lattice polyno- 
mials: 

Mx)= ^h{A) f\x, (xG [0,1]"), (3) 

AC[n.] ieA 

where rrih : 2^"' ^ M is the Mobius transform of Vh, defined as 

mh{A) := 5:(-l)l^l-l^lt;,(i?). 

BCA 

Indeed, expression ([3]) reduces to a linear function on each canonical simplex 
and agrees with h^ls) at 1b for each B <Z [n]. 
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Eq. ([2]) thus provides a concise expression for linear combinations of lattice 
polynomials. We will use it in the next section to calculate their distribution 
functions and their moments. 

Remark 1 ^45 we have already mentioned, the class of linear combinations 
of lattice polynomials covers three interesting particular cases, namely: lattice 
polynomials, linear combinations of order statistics, and weighted sums. These 
are characterized as follows. Let h : [0, 1]" — > M 6e a linear combination of 
lattice polynomials. 

(1) The function h reduces to a lattice polynomial if and only if the set func- 
tion Vh is monotone, {0, l}-valued, and such that Vh{[n]) = 1. 

(2) As the order statistics are exactly the symmetric lattice polynomials (see 
[n^), the function h reduces to a linear combination of order statistics 
if and only if the set function Vh is cardinality- based, that is, such that 
Vh{A) = VhiA') whenever \ A\ = \ A'\. 

(3) The function h reduces to a weighted sum if and only if the set function 
Vh is additive, that is, Vh{A) = I]ie/i ^/^({0)■ 



3 Distribution functions and moments 



Before yielding the main results, let us recall some basic material related to 
divided differences. See for instance [8]|ll|22j for further details. 

Consider the plus (resp. minus) truncated power function (resp. x"), de- 
fined to be if x > (resp. x < 0) and zero otherwise. Let A^^^ be the set of 
77,-1 times differentiable one-place functions g such that absolutely 
continuous. The nth. divided difference of a function g E A^"^^ is the symmetric 
function of n + 1 arguments defined inductively hj A[g : Oq] := g{aQ) and 

A[g : ai, . . . ,an] - A[g : ao, . . . ,an-i] , 
, II ao T ^ni 

A[(?:ao,...,a„] := { q ~ 

— : ao, . . . ,a„_i], if Oq = a„. 



The Peano representation of the divided differences, which can be obtained 
by a Taylor expansion of ^f, is given by 

A[g : ao, . . . , a„] = ^ / g^^\t) M{t | ao, . . . , a„) dt, (4) 



where M{t \ oq, . . . , a„) is the B-spline of order n, with knots {ao, . . . , a„}, 
defined as 

M{t I ao,...,a„) := n A[{- - t)l'^ :ao,...,a„]. (5) 
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We also recall the Hermite-Genocchi formula: For any function g G A^^\ we 
have 



/" r 

A[g : ao, . . . , a„] = / ^f^") Oq + E(«i " 



dx, (6) 



where Sid is the simplex defined in (|T]) when a is the identity permutation. 

For distinct arguments ao, . . . , a„, we also have the following formula, which 
can be verified by induction, 

A[g:a,,...,a.] = ± '[f (7) 



Now, consider a random vector X uniformly distributed on [0, 1]*^ and set 
Yh := /i(X), where the function h : [0,1]" — > M is a linear combination of 
lattice polynomials as given in formula ([2]). We then have the following result. 

Theorem 2 For any function g E A^"\ we have 

E[^?("Hr,)]= Y: A[g:K,...,hl]. (8) 



Proof. Using ([2]), we simply have 



E[g^-\Yh)]= f ^7(")[Mx)]dx 

"'[0,1]" 

^^(Tr 'J Srr L „■ 1 



dx. 



Finally, after an elementary change of variables, we conclude by the Hermite- 
Genocchi formula (El). □ 



Theorem [2] provides the expectation E[5r*^"^(Y/i)] in terms of the divided dif- 
ferences of g with arguments Zip, . . . , /i^ (o" G An explicit formula can be 
obtained by ([7]) whenever the arguments are distinct for every a G (5„. 

Clearly, the special cases 

^?(^) = r4^l r4^l - and ^ (9) 

[n + ry. [n + ry. t" 

give, respectively, the raw moments, the central moments, and the moment- 
generating function of Yh- As far as the raw moments are concerned, we have 
the following result. 



5 



Proposition 3 For any integer r ^ 1, setting Aq := [n], we have, 

[ r ) ^i^W -1 [ \M ) 

ArCAj,_i 



Proof. Let r ^ 1. It can be shown [3] that 

A[(.r-ao,...,a„]= Y: ■■■<"= E 

r-oH \-rn=r 

Hence, from ([S]) and (jH]) it follows that 



Era = 7;^ E E/^n---/^t 

= E E 

where Ain is the set of the n! maximal chains of the lattice (2["1,C), and 
where, for any m G A^„, mj is the unique element of m of cardinality i. 

For any i?i C ■ ■ ■ C 1?^ C [n], let ^Vl^i' - '^'' denote the subset of maximal 
chains of (2^"', C) containing Bi, . . . , Br- It is then easy to see that, for any 
fixed ^ ii ^ ■ ■ ■ ^ ir ^ n, the following identity holds: 

U A^fi' = Mn 

BiC---CBrC[n] 

\Bi\=ii,...,\Br\=ir 

and the union is disjoint. Therefore, we have 



Era = 7-^ E E E v,iB,)...v,iB..) 

' J- 0^ii^-<V^n BiC-CBrC[n] m&M^^'-''^'' 

Bl|=n,...,|i?r|=ir 



i=l 



{n + r)\ ^^c---cBrC[n,] 
where 

|A^f^'-'^1 = \B,\\ (\B,\ - \B,\)\ (IS3I - IS2I)! •■■(« - 
Finally, we get the result by setting Ai := Br+i-i for alH = 1, . . . , r. □ 
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Proposition [3] provides an explicit expression for the rth raw moment of Yh as 
a sum of (r + 1)" terms. For instance, the first two moments are 



"'^ ^ AC[n] ' ■ ■ ' 



(a,) 



, l^il 



We now yield a formula for the distribution function Fh{y) := Pr[Fft ^ y] 

of y/.. 

Theorem 4 There holds 

Fh{y) = i-^ E A[(.-y)!^:/.^,..., (10) 



Proof. We have 

Fh{y) = 1 - Pr[/i(X) > = 1 - E[(n - y)l 
Then, using ([8]) with 

leads to the result. □ 



It follows from ffTOl) that the distribution function of Yh is absolutely continu- 
ous and hence its probability density function is simply given by 

My) = E Mi-- y)l-' --K,---, K] (n) 

or, using the B-spline notation (jS]), 

fhiy) = -, E Miy\h-,,...,hl). 

fJGSn 

Remark 5 (i) It is easy to see that [T^) can be rewritten by means of the 
minus truncated power function as 

Fh{y) = ^ E ^{■-y)l:K,---,h:]- 

o-e&n 
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(ii) When the arguments h^, . . . ,h'^ are distinct for every a G ©„, then com- 
bining ^ with ( fl^) immediately yields the following explicit expression 



or, using the minus truncated power function, 



(Hi) The knowledge of fh{y) immediately gives an alternative proof of 
Indeed, using Peano's representation 0), we simply have 



E[^7("Hn)]= / g'^-\y) fniy) dy 
Jr 

nl ~z Jr 
= Y: A[g:h-„...,h:]. 

(iv) The case of linear combinations of order statistics is of particular interest. 
In this case, each h" is independent of a (see RemarklJ\), so that we can 
write hi := . The main formulas then reduce to (see for instance fJ^ 
and Ifj 

E[^7(")(n)]=n!A[^:/zo,...,M, 

Ff,{y) = A[{.-yr_:ho,...,hr,], 
fh{y) = M{y I /io,...,/i„). 

We also note that the Hermite- Genocchi formula ^ provides nice geo- 
metric interpretations of Fh{y) and fh{y) in terms of volumes of slices 
and sections of canonical simplices (see also /^i/ and 'il2j). 

Both functions Fhiij) and fh{y) require the computation of divided differences 
of truncated power functions. On this issue, we recall a recurrence equation, 
due to de Boor [9j and rediscovered independently by Varsi [21] (see also [3]), 
which allows to compute A[(- — y)""^ : ao, . . . , a„] in Oiji^) time. 

Rename as hi, . . . ,br the elements Cj such that Oi < y and clS C\ , . . . , Cg the 
elements such that ^ y so that r + s = n + 1. Then, the unique solution 
of the recurrence equation 

(q - y)ak-i,i + iy- bk)ak,i-i . , . ^ 

ak^i = [k^r, I ^ s), (12) 

ci — Ok 
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with initial values ai^i = (ci — 61) ^ and ao,/ = «fc,o = for all l,k ^ 2, is 
given by 

ak,i := A[(- - y)l+^'^ : h, . . . ,bk,Ci, . . . ,Ci] {k + l^2). 

In order to compute A[(- — y)"^^ : ao, . . . ,an] = ar,s, it suffices therefore to 
compute the sequence ak,i for k + 1 ^ 2, k ^ r, I ^ s, hj means of two nested 
loops, one on k, the other on /. 

We can compute A[{- — y)^ : ao, . . . , an] similarly. Indeed, the same recurrence 
equation applied to the initial values ao,; = for all I ^ 1 and akfi = 1 for all 
k ^ 1, produces the solution 

ak,i := A[(. - yt^'-' : 61, . . . , 6^, Ci, . . . , q] (k + l^l). 

See for instance [1] and [21] for further details. 



4 Application to aggregation theory 

As we have already mentioned, the concept of linear combination of lattice 
polynomials, when it is nondecreasing in each variable, is known in aggrega- 
tion theory as the discrete Choquet integral, which is extensively used in non- 
additive expected utility theory, cooperative game theory, complexity analysis, 
measure theory, etc. (see [13] for an overview.) For instance, when a discrete 
Choquet integral is used as an aggregation tool in a given decision making 
problem, it is then very informative for the decision maker to know its dis- 
tribution. In that context, the most natural a priori density on [0, 1]" is the 
uniform one, which makes the results derived here of particular interest. 

Example 6 Let h : [0, 1]'^ M 6e a discrete Choquet integral defined by 
Vh{{l}) = 0.1, vh{{2}) = 0.6, v^m) = vh{{1, 2}) = 3}) = Vh{{2, 3}) = 

0.9; and Vh{{l,2,3}) = 1. According to it can be written as 

/i(x) = 0.1 xi + 0.6 X2 + 0.9 X3 

+ 0.2{xi A X2) - 0.l{xi A X3) - 0.6(x2 A X3) 
— 0.1(zi A 2:2 A X3). 

Its density, which can be computed through 07]) and the recurrence equa- 
tion ( fl^) . is represented in FigureUl by the solid line. The dotted line represents 
the density estimated by the kernel method from 10 000 randomly generated 
realizations. The typical value and standard deviation can also be calculated 
through the raw moments: we have 

E[Yfe] ^ 0.608 and ^^[Y^] - E[n]2 ^ 0.204. 
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Fig. 1. Density of a discrete Choquet integral (solid line). 
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